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Origin of ozone and NO x in the tropical troposphere: A 
photochemical analysis of aircraft observations over the 
South Atlantic basin 

D.J. Jacob , 1 B.G. Heikes , 2 S.-M. Fan , 1 ’ 3 J.A. Logan , 1 D.L. Mauzerall , 1 J.D. 
Bradshaw , 4 H.B. Singh , 5 G.L. Gregory , 6 R.W. Talbot , 7 D.R. Blake , 8 and G.W. 
Sachse , 6 

Abstract. The photochemistry of the troposphere over the South Atlantic basin is 
examined by modeling of aircraft observations up to 12 -km altitude taken during the 
TRACE A expedition in September-October 1992. A close balance is found in the 0 to 
12 -km column between photochemical production and loss of CL, with net production at 
high altitudes compensating for weak net loss at low altitudes. This balance implies that 
CL concentrations in the 0-12 km column can be explained solely by in situ 
photochemistry; influx from the stratosphere is negligible. Simulation of H : 0 : . 
CFLOOH, and CH 2 0 concentrations measured aboard the aircraft lends confidence in 
the computations of O 3 production and loss rates, although there appears to be a major 
gap in current understanding of CH 2 0 chemistry in the marine boundary layer. The 
primary sources of NO x over the South Atlantic Basin appear to be continental (biomass 
burning, lightning, soils). There is evidence that NO x throughout the 0 to 12-km column 
is recycled from its oxidation products rather than directly transported from its primary 
sources. There is also evidence for rapid conversion of HNCL to NO x in the upper 
troposphere by a mechanism not included in current models. A general representation of 
the 0 3 budget in the tropical troposphere is proposed that couples the large-scale Walker 
circulation and in situ photochemistry. Deep convection in the rising branches ol the 
Walker circulation injects NO x from combustion, soils, and lightning to the upper 
troposphere, leading to O 3 production; eventually, the air subsides and net CL loss takes 
place in the lower troposphere, closing the O 3 cycle. This scheme implies a great 
sensitivity of the oxidizing power of the atmosphere to NO x emissions in the tropics. 


1. Introduction 

Ozone in the tropical troposphere plays a key role in deter- 
mining the global oxidizing power of the atmosphere. Most of 
the oxidation of long-lived gases by OH takes place in the tropics, 
where high UV and humidity promote the formation of OH from 
photolysis of 0 3 [Logan et al , 1981; Thompson , 1992]. Ozone is 
produced within the troposphere by photochemical oxidation of 
hydrocarbons and CO in the presence of nitrogen oxides (NO v = 
NO + N0 2 ) and is also transported down from the stratosphere. 
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There is continuing controversy as to the relative importance of 
tropospheric production versus stratospheric influx in controlling 
Ox concentrations in the troposphere, including in the tropics [Lin 
et al., 1980; Levy et al., 1985; Logan and Kirchhoff, 1986; Krish - 
naniurti et al.. 1993; Dibb et ai, 1996; Graustein and Turekian , 
1996; Browed et a /., this issue; Loring et al. , this issue]. There is 
also considerable uncertainty regarding the origin of NO x , the 
limiting precursor for CL production in the troposphere [Levy and 
Max ini . 1989; Penner et al . , 1991; Kasibhatla et al. . 1991, 1 993 ]. 

We examine here the origin of CL and NO v in the tropical tro- 
posphere using observations made during the TRACE A aircraft 
expedition in September-October 1992 [Fishman et al., this 
issue]. This expedition provided a detailed survey of and its 
precursors over the tropical South Atlantic and neighboring con- 
tinents up to 12-km altitude. Its primary objective was to identify 
the cause of the large-scale O* enhancement observed in the 
region during the austral spring dry season [Fishman et al.. 1 990 1. 
Results from TRACE A confirmed earlier suggestions that this 
enhancement is caused by seasonal biomass burning emissions of 
NO x , hydrocarbons, and CO from southern Africa and South 
America [Logan and Kirchhoff. 1986; Fishman et al.. 1990; 
Anderson et al., 1993; Andreae et al., 1994; Gregory et al., this 
issue; Kirchhoff et al., this issue; Olson et al.. this issue; Thomp- 
son et ai, this issue]. 

Our goal in the present paper is to explain the total tropos- 
pheric 0 ; column observed in TRACE A. not just the biomass 
burning enhancement, and to draw general inferences regarding 
[he budget of Oj in the tropical troposphere. As pointed out by 
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Thompson et ai [this issue], the tropospheric 0 } column over the 
South Atlantic during the dry season maximum is only 50% 
higher than the wet season background when biomass burning 
influence is minimal. This background, over the South Atlan- 
tic and elsewhere in the tropical troposphere, plays a critical role 
in determining the oxidizing power of the atmosphere. 

Our analysis is based on a combination of photochemical 
modeling and statistical interpretation of the TRACE A data. The 
modeling approach is described in section 2 and is evaluated in 
sections 3 and 4 by simulation of H 2 0 2 , CH^OOH, and CH 2 0 
concentrations measured aboard the aircraft. Regional budgets 
for 0 3 and NO x over the South Atlantic basin are derived in sec- 
tions 5 and 6, respectively. A general representation of the O, 
budget in the tropical troposphere is proposed in section 7. Con- 
clusions are in section 8. 

2. Modeling Approach 

The TRACE A expedition used a DC-8 aircraft with a ceiling 
of 12-km altitude. The southern hemisphere flight tracks (Figure 
1) covered the tropical South Atlantic, eastern Brazil, and south- 
ern Africa; we refer to this region as the 'South Atlantic basin. 
Details on the flights and references for the observations are 
given by Fishman et ai [this issue]. Measurements aboard the 
aircraft included concentrations of 0 3 , CO, NO, peroxyacetylni- 
trate (PAN), peroxypropionylnitrate (PPN), HN0 3 , speciated 
hydrocarbons, acetone, H 2 0 2 , CH3OOH, CH 2 0, and organic 
acids. Meteorological variables, UV fluxes, and aerosol number 



Figure 1. TRACE A flight tracks (September 23 - October 25, 
1992). Flight numbers are indicated. The region covered by the 
southern hemisphere flight tracks is referred to in the paper as the 
South Atlantic basin. 


size distributions were also measured. Table 1 gives median 
values for the various measurements at different altitudes. 

We base our analyses on merged time series of the aircraft 
data where all measurements have been averaged over a common 
time interval. The choice of averaging interval must balance the 
need for high temporal resolution with the danger of misinterpret- 


Table la. Median Measured Values Over the South Atlantic Basin During TRACE A 




Altitude 



0-4 km 

4-8 km 

8-12 km 

Temperature, K 

284 

258 

230 

Dew point, K 
UV flux. a mW cm " 2 

275 

244 

218 

zenith 

4.7 

5,3 

5.8 

nadir 

0.75 

1.5 

2.0 

Oi. ppbv 

50 

69 

74 

CO, ppbv 

112 

103 

93 

CH 4 . ppbv 

1695 

1707 

1715 

NO, pptv 

25 

29 

120 

PAN. pptv 

175 

294 

223 

PPN, pptv 

4 

3 

2 

HNOv pptv 

390 

130 

59 

Ethane, pptv 

590 

630 

670 

Propane, pptv 

48 

46 

58 

C 4.5 alkanes, pptv 

9 

7 

8 

Ethene, pptv 

21 

8 

9 

Propene. pptv 

6 

4 

4 

Acetylene, pptv 

130 

130 

1 10 

Benzene, pptv 

32 

21 

15 

Toluene, pptv 

7 

2 

<2 

Acetone, pptv 

750 

593 

650 

Ethanol, pptv 

60 

<20 

<20 

HCOOH, pptv 

1770 

1000 

560 

CH^COOH, pptv 

3030 

2550 

1900 

H 2 0 2 , pptv 

2800 

760 

150 

CH^OOH. pptv 

1200 

320 

60 

CH 2 0, pptv 
Aerosol number, cm ' 1 

200 

<40 

<40 

0.12 to 0 . 2 -gm radius 

181 

35 

16 

0.2 to 0.5-gm radius 

310 

28 

8 

0.5 to 2.0-pm radius 

0.2 

< 0.1 

0.3 


Statistics are for the ensemble of southern hemisphere data collected in TRACE A. The median 
concentrations of >C 5 alkanes, >Cy alkenes, xylenes, isoprene, and terpenes were below the detec- 
tion limit of 2 pptv. Concentrations of N0 2 were measured in TRACE A. but the measurement is 
considered unreliable [ Crawford et ai, 1996] and is not used in the paper. 

a Measured with an Eppley radiometer (Eppley, Inc., Newport, Rhode Island). 
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Table lb. Median Model Values Over the South Atlantic Basin During TRACE A 




Altitude 


0-4 km 

4-8 km 

X. 12 km 

OH. 10 s cm' 1 

27 

16 

12 

HO->, I0 6 cm' 1 

257 

96 

22 

CH,0 2 . I0 6 cm' 1 

103 

34 

3 

NO,=NO+NO->, pptv 

73 

57 

145 

HNOj, pptv 

1 

10 

54 


Twenty-four hour average model values computed as described in section 2. The concentration of 
NO* is calculated in the model using as constraint the NO measurement for the time of day of obser- 
vation. 


ing coarsc-resolution data. We choose to average the data over 
the sampling time of the critical variable with the coarsest time 
resolution. Because this critical variable depends on the problem 
at hand, we will use over the course of this paper four different 
merged time series matched to the sampling times of peroxides, 
CH 2 0, NO, and HNO, (Table 2). Ozone columns are included in 
the merged time series using daily observations with l°xl.25° 
resolution from the total ozone mapping spectrometer (TOMS). 
The median O, column over the south tropical Atlantic during 
TRACE A was 281 Dobson units (DU) (1 DU = 2.687x1 0 lfl 
molecules cm' 2 ). 

A major objective of this paper is to derive a chemical budget 
for O, in the tropospheric column over the South Atlantic basin 
(section 5). Our approach is to calculate 24-hour average O, pro- 
duction and loss rates for individual points in the merged time 
series of aircraft data and to use statistics from the ensemble of 
points to infer regionally averaged rates. Rapid chemical cycling 
takes place between O, and other members of the odd-oxygen 
family (O x = O, + O + N0 2 + HN0 4 + (2x)NO, + (3x)N 2 O s + 
pemxyacetylnitrates), hence the budget analysis must be based on 
O y rather than on O v Considering that O, accounts typically for 
over 99% of O x , the budgets of O, and O x can be viewed as 
equivalent. Production of O x is mainly by reactions (1) and (2): 

H0 2 + N0->0H+N0 2 (1) 

R0 2 + NO -> RO + N0 2 (2) 

where R0 2 is any organic peroxy radical (CH,0 2 is the most 
important). Loss of O x is principally by reactions ,(3)-(5): 

O, + hv -> 0 2 + 0( 1 D) (3a) 

0(‘ D) + H 2 0 —> 20H (3b) 


H0 2 + O, — > 20H (4) 

OH + O, — > H0 2 4- 0 2 (5) 

Calculation of O x production and loss rates requires 
knowledge of the concentrations of H0 X radicals (OH, H0 2 , 
R0 2 ) along the flight tracks. These species were not measured 
from the aircraft. We estimate their concentrations for each point 
in the merged time series by using a zero-dimensional photo- 
chemical model constrained with the ensemble of measurements 
for that point (Table 3). The model is described in Appendix A. 
It calculates the diel steady state concentrations of H0 X radicals, 
short-lived NO y species (NO, N0 2 , NO,, N 2 0 5 , HN0 2 , HN0 4 , 
organic nitrates), and other photochemical intermediates (perox- 
ides, carbonyls). Production and loss of species are by chemical 
reactions only (loss by deposition is not considered). "Diel 
steady state" is defined by reproducibility of concentrations over 
a 24-hour solar cycle, as obtained in a time-dependent calculation 
with periodic boundary conditions 

C,(0 = C,(r+77 (6) 

where C, is the concentration of species 1 , t is time, and T - 24 
hours. In this calculation, all input variables in Table 3 except 
NO are assumed constant over the diel cycle. The concentration 
ot NO must be calculated as a time-dependent quantity because 
of the strong diel variation driven by photochemical steady state 
with N0 2 . The few available observations in the remote tropo- 
sphere indicate little diel variation of N0 X [Carroll et a/., 1992], 
but some variation might be expected in the upper troposphere 
due to cycling with HN0 4 . We assume in the model that the con- 
centration of N0 t (NO + N0 2 + NO, + (2x)N 2 0 5 + HN0 2 + 
HN0 4 ) is constant over the diel cycle and calculate its value 
within the model to match the observed NO concentration at the 
time of day of the measurement. 


Table 2. Merged Time Series of TRACE A Data 


Matched Species 

Time 

Resolution, 

min 

Data 

Number of Points 
Model 

Application 

Peroxides 

4-5 

1254 

223 

a 

ch 2 o 

1 

1496 

61 (926) 

b 

NO 

3 

1781 

214(1355) 

c 

HNO, 

2-15 

721 

201 

d 


Each data point represents one measurement time interval for the matched species. Model points 
are the subset of data points for which photochemical model quantities can be computed; the compu- 
tation requires that data for all variables in Table 3 be available for at least a fraction of the time 
interval. We also give in parentheses the number of model points obtained in the fixed -hydrocarbon 
approximation, which relaxes the constraint on availability of data for hydrocarbons and acetone (see 
section 2). 

:i Comparison of simulated and observed peroxide concentrations (section 3). 
h Comparison of simulated and observed CH 2 0 concentrations (section 4). 
c O, production and loss rates and N0 2 concentrations (sections 5 and 6). 

J Chemical budget of N0 X (section 6). 
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Table 3. Input Variables for the Photochemical Model Calculation 


Variable Type 

Variable 

Meteorological 

pressure, temperature, humidity 

Radiation 

solar declination, latitude, 0; column, J 
cloud correction factor CF h 

Chemical 

0;,CO, NO c .PAN d ,HNOi 

methane 0 , ethane, propane, C 4 _5 alkanes, Q ,_7 

alkanes, ethene, propene, C 4 _* alkenes, acetylene, benzene, 

toluene, xylenes, isoprene, acetone, ethanol, formic acid'' , acetic acid . 

aerosol number concentrations in three size classes'- 1 


a The total O; column is distributed vertically using the mean O; profiles measured by ozo- 
nesondes at Ascension Island and southern Africa during TRACE A [Thompson et «/., 1995|. 

h The cloud correction factor (CF) is calculated from the UV flux observations as described in 


c The model calculation assumes a constant concentration of NO t (NO, - NO + NO : + NO; + 
(2 x)N 2 O s + HN0 2 + HN0 4 ) over the diet cycle and adjusts the concentration of NO, to match the 

observed NO concentration at the time of day of observation. . 

d Data for HNO; and PAN are required only for the NO x budget calculation ot section 6. For 
other applications, if HNO; data are unavailable then a low value is assumed: and if PAN data are 
unavailable, then PAN concentrations are assumed from chemical steady state. These assumptions 
greatly increase the number of data points at which model calculations can be done without inducing 
significant error in the relevant model output quantities. 

c If data are not available for these variables, default values are used ot 1.7 ppmv tor CH 4 and 
zero for formic acid, acetic acid, and aerosol number concentrations. 

1 See Table la for size class information. 


Concentrations of H 2 0 2 , CH;OOH, and CH 2 0 are calculated 
in the model, rather than constrained from observations, and com- 
parison of simulated versus observed values serves to test the 
model (sections 3 and 4). These three species are sensitive to dif- 
ferent aspects of HO x chemistry [Heikes, 1992; Liu et al . , 1992]. 
Hydrogen peroxide is particularly valuable as a test ot O; produc- 
tion in the model because its source, 

H0 2 4- H0 2 -> H 2 0 2 + 0 2 (7) 

varies quadratically with the concentration ot H0 2 , the principal 
peroxy radical responsible for O; production. At low altitudes 
where H 2 0 2 is the principal sink of HO x , the concentration of 
H,0 2 reflects the rate of the 0( l D)+H 2 0 reaction which is the 
principal source of HO x and the principal sink of O x . 

We will also apply the model to study the chemical cycling of 
NO x over the South Atlantic basin in light of the concurrent 
measurements of NO, HNO;, and PAN made aboard the aircraft 
(section 6). Oxidation of NO x to HNO; and PAN takes place on 
a timescale of the order of 1 day: 


N0 2 + OH + M -> HNO, + M 

(8) 

N0 2 + 0, -> NO, + 0 2 

(9a) 

NO, + NO, + M -> N 2 0 5 + M 

(9b) 

aerosol. H,0 

N 2 O s - J —> 2HNO, 

(9c) 

NO, + CH,CO, + M — » PAN + M 

(10) 

Organic nitrates other than PAN are also formed but represent 
only minor sinks for NO x in the model. The species HNO; and 
PAN can act as reservoirs for NO x , decomposing to release NO x 
far from its point of origin: 

HNO, + hv -> N0 2 + OH 

(H) 

HNO, + OH -> NO, + H 2 0 

(12a) 


NO; + hv -» N0 2 + O 

(12b) 

PAN heal >NO, +CH,CO, 

(13) 

PAN + hv -» NO, + CH,CO, 

(14) 


We compute the 24-hour average rates of reactions (8)-( 14) in the 
model at individual points using the observed concentrations of 
NO, HNO; and PAN as constraints, as discussed above. The 
resulting mass balance between sources and sinks of NO x allows 
us to estimate the degree to which NO x concentrations are main- 
tained by chemical recycling from HNO; and PAN through the 
above set of reactions. In such a NO x budget calculation it is 
necessary to correct the rates of conversion between NO x and 
PAN to account for cycling within the (CH;CO;+PAN) family: 
the correction is described in Appendix B. 

Table lb gives the median concentrations of OH. H0 2 . 
CH ; Oi , NO x , and HN0 4 computed in the model for the southern 
hemisphere TRACE A data set. Diel steady state tor these 
species and for CH 2 0 is typically approached to within 10% after 
2 days of integration in the model. The peroxides, which have 
longer lifetimes, require typically 2-3 days at 2 km and 1 week at 
8-12 km for approach to steady state. The steady state assump- 
tion is prone to large errors in pollution plumes and in the con- 
tinental boundary layer, where concentrations ot the input vari- 
ables in Table 3 may change greatly over timescales of less than I 
day. It is more adequate in the free troposphere and in the marine 
boundary layer. 

To conduct the model calculation for a given point in the 
merged time series we require that data be available tor all vari- 
ables in Table 3. Because of lack of temporal overlap in some of 
the measurements (notably NO, hydrocarbons, acetone), model 
calculations can be done for only a small fraction of the data 
points (Table 2). To increase the number of points, we use in cer- 
tain applications a "fixed-hydrocarbon” version of the model in 
which we assume 630 parts per trillion by volume (pptv) C 2 H 6 
and 650 pptv acetone (mean values in TRACE A) and neglect the 
other nonmethane hydrocarbons. By relaxing the constraint on 
the availability of hydrocarbon and acetone data, the number of 
model points is greatly increased (Table 2). The errors incurred 
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by the fixed hydrocarbon approximation will be discussed in the 
context of each application. Data for HNO^ and PAN are critical 
only for the NO x mass balance analysis presented in section 6. 
For all other applications, when data for HNO^ are missing a low 
value is assumed; and when data for PAN are missing the PAN 
concentration is assumed from chemical steady state. We verified 
that these assumptions have a negligible effect on HO x levels and 
other relevant model results. 

3. Peroxides 

Figure 2 compares simulated and observed concentrations of 
H 2 0 2 and CH^OOH for individual points in the merged time 
series matched to the peroxide sampling times. Model values 
were sampled from the diel steady state cycle at the time of day 
of observation. Two outliers representing fresh fire plumes were 
excluded from the plots and from the statistics. Also shown on 
each plot are the least squares correlation coefficient r between 
model and observations, the slope S of the linear best fit obtained 
by the reduced-major-axis method [Hirsch and Gilroy , 1984], and 
the median ratio R of simulated-to-observed concentrations. S 
and R offer different measures of model bias. 

We find that the model captures 78% of the observed variance 
of H 2 0 2 concentrations. This variance is determined strongly by 
water vapor (r = 0.51 between observed H 2 0 2 and H 2 0 mixing 
ratios). Median simulated-to-observed concentration ratios are 
0.99 at 0-4 km, 0.87 at 4-8 km, and 0.57 at 8-12 km. Although 
H 2 0 2 deposition is not included in the model, its effect should be 
largely limited to the mixed layer and to convective outflows 
because the lifetime of H 2 0 2 against photochemical loss is rela- 
tively short (1.5 days below 2 km). The underestimate of H 2 0 2 




OBSERVED, ppbv 

Figure 2. Simulated versus observed H 2 0 2 and CH^OOH concen- 
trations in TRACE A. Also shown are the 1:1 line, the least 
squares correlation coefficient (r ), the slope of the regression line 
(5), and the median ratio of simulated to observed concentrations 
(/?). Observed CH^OOH concentrations were below the detection 
limit of 15 pptv for 15% of the points. 



0.01 0.1 1 
OBSERVED, ppbv 

Figure 3. Simulated versus observed CH 2 0 concentrations in 
TRACE A. Also shown are the 1:1 line, the least squares correla- 
tion coefficient (r), the slope of the regression line (5), and the 
median ratio of simulated to observed concentrations (R). 
Observed concentrations were below the detection limit of 30-50 
pptv for 56% of the points. 


in the upper troposphere could reflect a missing HO x source in 
the model but is also within the low-temperature uncertainties of 
the rate constants in the mechanism. 

For CH^OOH we find that the model captures 74% of the 
observed variance but is generally too low, particularly above 8 
km. Median simulated-to-observed ratios are 0.59 at 0-4 km. 
0.79 at 4-8 km. and 0.28 at 8-12 km. The underestimate is not 
outside the range of uncertainty of rate constants in the mechan- 
ism [Thompson and Stewart , 1991]. The estimated uncertainty on 
the rate constant for the reaction CH^+HOt, which is the only 
CH^OOH source in the model, is a factor of 2 at 298 K and a fac- 
tor of 3 at 230 K [DeMore et al . , 1994]. Although the underesti- 
mate of CHiOOH could suggest the presence of major hydrocar- 
bon precursors missing from the model, we find no parallel 
underestimate of CH 2 0 that would support this hypothesis (see 
section 4). 

The successful simulation of the peroxides below 8-km alti- 
tude lends confidence in the 0 3 production and loss rates com- 
puted by the model. At 8 to 12-km altitude the ~ 40% underesti- 
mate of H 2 0 2 could imply an - 20% underestimate of HO : . 
Considering that H0 2 accounts for -90% of total peroxy radicals 
in that altitude range (Table lb), we may expect the 0 3 produc- 
tion and loss rates above 8 km to be underestimated by about 
20 %. 

Singh et al. [1995] pointed out the importance of acetone pho- 
tolysis (included in our model) as a primary source of HO s in the 
upper troposphere. We conducted a sensitivity calculation 
without acetone for the point representative of median conditions 
at 8-12 km (Table la). The H 2 0 2 , CH 3 OOH, H0 2 , and OH con- 
centrations computed for that point decreased by 27% , 36%, 
14%', and 16%, respectively, relative to the simulation including 
acetone. Neglecting acetone as a source of HO x makes the model 
underestimate of peroxide concentrations in the upper tropo- 
sphere significantly worse. 

4. Formaldehyde 

Figure 3 compares simulated and observed CH 2 0 concentra- 
tions for individual points in the merged time series matched to 
the CH 2 0 sampling times. Again, model values were sampled 
from the diel steady state cycle at the time of day of observation. 
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The model captures 62% of the observed variance. Model values 
tend to be too high in the lower troposphere (0-4 km) and too low 
in the middle troposphere (4-8 km). Most observations above 8 
km are below the detection limit of 40 pptv. 

Our ability to interpret model results for CH 2 0 is evidently 
hampered by the small number of model points. To address this 
problem, we turned to the fixed-hydrocarbon version of the model 
which provides 926 model points for comparison with observa- 
tions. The fixed-hydrocarbon approximation gives only a lower 
limit for CH 2 0, because it does not allow for the possibility of 
high hydrocarbon concentrations, but it still allows investigation 
ofmodel overestimates. We find that the worst overestimate is in 
the marine boundary layer over the tropical South Atlantic (mean 
observed conditions 26 parts per billion by volume (ppbv) CL, 66 
ppbv CO, 3 pptv NO, 0.015 vol/vol H 2 0), where the mean 
observed concentration of CH.O is 0.1 1 ppbv (range <0.04-0.25 
ppbv) while the corresponding model mean is 0.41 ppbv (range 
0.28-0.55 ppbv). The mean observed H 2 0 2 and CH 3 OOH con- 
centrations in that region are 2.3 ppbv and 1.1 ppbv, respectively, 
while the corresponding model means are 3.1 ppbv and 1.7 ppbv. 
The model shows some overestimate for the peroxides, which in 
the case of H 2 0 2 could be due to deposition [Heikes et ai , this 
issue], but these discrepancies are small compared to CH 2 0. We 
repeated the model calculation for the marine boundary layer 
points using observed concentrations of H 2 0 2 and CH 3 OOH as 
constraints, and assuming a deposition rate constant ot 4xKT* s 1 
for CH.O as representative of a deposition velocity of 0.4 cm s 1 
over the ocean [Thompson and Zafiriou , 1983]. The mean CH 2 0 
concentration in this calculation was 0.29 ppbv (range 0.21-0.48 
ppbv), still 3 times higher than observed. The root of the 
discrepancy must therefore lie in the simulation of CH.O chemis- 
try. 

Our model results for CH 2 0 can be viewed as typical ot the 
current generation of photochemical models. We participated 
recently in an intercomparison of 20 models organized by the 
Intergovernmental Panel on Climate Change [Prather et at ., 
1995]. In a test simulation representative of the remote marine 
boundary layer we obtained a CH 2 0 concentration of 0.21 ppbv, 
while other models obtained values ranging from 0.13 to 0.33 
ppbv. The models at the low range had unusually high CH.O 
photolysis frequencies. 

We compared the marine boundary layer measurements of 
CH.O in TRACE A to prior shipboard measurements made over 
the tropical Atlantic in the same season. Lowe and Schmidt 
[1983] measured a mean concentration of 0.20 ppbv (range 0.10- 
0.34 ppbv) at 30°N-30°S latitude with no significant latitudinal 
dependence. These values are higher than observed in TRACE A 
but still lower than in our model. Harris et ai [1992] used tun- 
able diode laser spectroscopy to measure CH.O concentrations in 
the northern tropics during the Polarstern cruise and found a 
mean concentration of 0.37 ppbv at 10°-20°N (they made no 
measurements farther south). During the same cruise, Carlier et 
cd . [1991] reported mean CH.O concentrations of 0.4 ppbv and 
0.3 ppbv over the north and south tropical Atlantic, respectively. 
We evaluated our model against the Harris et ai [1992] data by 
taking as constraints the mean Polarstern observations at 15°N: 
30 ppbv 0 3 , 90 ppbv CO, 1.3 ppbv H 2 0 2 , 15 pptv NO., 0.017 
mol/mol H O, temperature of 293 K, and actinometer values of 
J N(h and j' () i D [Platt et ai , 1992; Jacob and Klockow , 1992; 
Bra uers and Hofzumahaus, 1992; Hofzumahaus et ai , 1992]. We 
obtained a model CH.O concentration of 0.42 ppbv, in good 
agreement with the observations. 

Other published measurements of CH.O concentrations in 
clean tropical marine boundary layers include 0.4- 1.1 ppbv over 
the Pacific [Zafiriou et ai , 1980; Arlander et ai , 1990] and 0.05- 
0.4 ppbv over the Indian Ocean [Arlander et ai , 1990]. It is not 
clear why there is such variability in the observations and why 
models are unable to reproduce the lower range of values. One 
possibility might be the complexation of CH.O with organics in 


the marine aerosol. To our knowledge, the only concurrent gas 
phase and aerosol measurements of CH 2 0 in marine air are those 
of Klippel and Warneck [1980], These authors found evidence 
that CH 2 0 in the aerosol is predominantly present as CH.O com- 
plexes, but they also found that the aerosol contributes only a 
small fraction of total atmospheric CH.O. 

In summary, our results point to a major gap in current under- 
standing of CH.O concentrations in the marine boundary layer. 
Either some of the CH 2 0 measurements are flawed or a large 
CH.O sink is missing from the models under some conditions. 
Liu et ai [1992] previously noted that CH 2 0 concentrations 
measured in free troposphere air at Mauna Loa, Hawaii (3.4-km 
altitude) were a factor of 3 lower than calculated from a standard 
photochemical model and attributed the discrepancy to a model 
overestimate of OH. However, later measurements by Zhou et ai 
[1996] at Mauna Loa indicated that the discrepancy is driven by a 
model overestimate of CH 3 OOH and that CH.O concentrations 
calculated in a model constrained with observed CH 3 OOH are in 
good agreement with observations. A global three-dimensional 
photochemical model study by Brasseur et ai [1996] shows, in 
fact, good simulation of both CH 3 OOH and CH 2 0 at Mauna Loa. 
We cannot draw a parallel between Mauna Loa and the marine 
boundary layer in TRACE A, because the discrepancy between 
model and observations for CH 2 0 in TRACE A is far greater 
than for CH 3 OOH; also, as discussed above, there is no evidence 
in TRACE A for a model overestimate of CH.O in the free tropo- 
sphere. We searched the TRACE A data set for conditions that 
would approximate the Mauna Loa environment described by Liu 
et ai [1992] and Zhou et ai [1996], but all TRACE A measure- 
ments at that altitude had much lower temperatures and humidi- 
ties than at Mauna Loa (temperature and humidity are important 
variables controlling CH.O concentrations in the model). 

5. Ozone 

Plates la and lb show the horizontal distributions of the 24- 
hour average gross production rate of odd oxygen, P ()k , and the 
net production rate, (P-L) 0x , at different altitudes in TRACE A. 
The rates were computed for the merged time series matched to 
the NO sampling time, using the fixed-hydrocarbon approxima- 
tion in order to maximize the number of model points. The errors 
on P a and L ()x caused by the fixed-hydrocarbon approximation 
are small; for those points where the full model (including hydro- 
carbons) can be applied, the median ratio of P 0x calculated in the 
full versus fixed-hydrocarbon versions of the model is 
1.02±0.09(1 g) when P Qx is less than 20 ppbv day" 1 . This regime 
represents 92% of all observations in TRACE A and all observa- 
tions outside the continental boundary layer. Plate lb indicates a 
broad region of net photochemical loss exceeding 4 ppbv day -1 at 
1.5 to 4-km altitude off the southwestern coast of Africa. This 
rapid loss is made possible by the albedo of persistent marine 
stratus decks. Gross photochemical production in the region is 
also high. 

We see from Plate 1 that the gross production of 0 3 decreases 
on average with altitude while the net production increases with 
altitude. Median vertical profiles of P Qx and L a% for the southern 
hemisphere TRACE A data are plotted in Figure 4, and a sum- 
mary 0 3 budget for the 0 to 12-km column is given in Table 4. 
Both production and loss decrease with altitude, due to decreas- 
ing humidity, but production decreases slower than loss above 6 
km because of high concentrations of NO. Above 6 km, photo- 
chemistry switches from a net sink to a net source of 0 3 . Similar 
transitions have been reported in previous model studies of the 
remote troposphere [Liu et ai , 1980; Thompson and Cicerone , 
1982; Davis et ai , 1996]. Figure 4 indicates net production of 0 3 
taking place below 2-km altitude, but this portion of the budget is 
highly uncertain because of the large variability of 0 3 production 
rates for individual points; there is a sharp contrast between the 
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24-Hour POx in ppb/day 

Plate la. Twenty-four-hour average gross production rates of odd 
oxygen (O x = O3 + O + NO2 + HNO4 + (2x)N03 + (3x)N20s + 
peroxyacetylnitrates) computed at 0 to 4-km, 4 to 8-km, and 8 to 
1 2-km altitude in TRACE A. 


continental boundary layer, where biomass burning influence can 
result in rapid O y production, and the marine boundary layer, 
where NO x levels are low and O y production is slow (Plate 1). 
The point modeling approach is also most prone to errors in the 
continental boundary layer, as discussed in section 2. 

We take the values in Figure 4 as our best estimates of the 
regionally averaged rates of O3 production and loss over the 
South Atlantic basin. The corresponding integrals of P ()x and L ()x 
for the 0 to 12-km column are both l.lxl 0 12 cm -2 s" 1 . These 
values are two orders of magnitude larger than the mean stratos- 
pheric influx of 0 3 in the southern hemisphere troposphere 
estimated to be l-3xl0 n> cm -2 s” ! [Gidel and Shapiro, 1980; Levy 
et al. , 1985]. Even allowing for a slightly higher cross- 
tropopause flux in spring f Holton , 1990] and for preferential sub- 
sidence of air over the south tropical Atlantic [Krishnamurti et 
al., 1993; Loring et al , this issue], it is clear that stratospheric 
influx is a negligible contributor to tropospheric 0-, in the region. 
In situ photochemistry dominates even in the upper troposphere; 
the 4-km integral of P 0x at 8 to 12-km altitude is l.lxlO 11 cm -2 
s" 1 . 

Deposition has only a small effect on the O3 budget in the 0 
to 12-km column. Heikes et al. [this issue] estimated a deposition 
flux to the Atlantic Ocean of 1.7xl0 1() cm -2 s' 1 during TRACE A. 
Fan et al. [1990] measured a 24-hour average O, deposition flux 


of 10xl0 M) cm" 2 s" 1 to the Amazon forest in the wet season; scal- 
ing this flux to the ratio of dry and wet season 0 } concentrations 
over the forest [Browell et ai f 1990] would imply a dry season 
deposition flux to the Amazon forest of 17xl0 ln cm" 2 s" 1 . 
Another term in the O3 budget is the time tendency, but this term 
is also small. Seasonal trends in the tropospheric O3 column at 
sites in the South Atlantic basin are of the order of 5 DU per 
month, or 5x1 0 10 molecules cm" 2 s" 1 [Olson et al., this issue]. 

The close balance between chemical production and chemical 
loss of O3 in the 0 to 12-km column implies that the total column 
amount of 0 } can be explained by in situ photochemistry, with 
net production above 6-km altitude balancing weak net loss 
below. To demonstrate this point, we conducted a one- 
dimensional photochemical steady state model calculation for O y 
in the 0 to 12-km column, fixing as input the median vertical 
profiles of O3 precursors and other conditions from Table la and 
adopting a zero-flux upper boundary condition at 12-km altitude. 
The calculation assumed an O3 deposition velocity of 0.05 cm 
s" 1 . a uniform vertical eddy diffusion coefficient K y = 3xl0 5 cm" 
and 24-hour average clear-sky photolysis rates. Results in 
Figure 5 show that it is possible to explain the O y concentrations 
and vertical gradients observed in the 0 to 12-km column in 
TRACE A without invoking any exogenous O y source. 

The median lifetime of O x against chemical loss over the 
South Atlantic basin is 6 days at 0-4 km, 25 days at 4-8 km, and 



24-Hour (P-L)Ox In ppb/day 

Plate lb. Twenty-four-hour average net O x production rates com- 
puted at 0 to 4-km, 4 to 8-km, and 8 to 12-km altitude in TRACE 
A. 
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Figure 4. (Left) Median 24-hour average production and loss rates of O, in the southern hemisphere TRACE A data as 
a function of altitude. (Right) Median concentrations of 0 3 and NO. 


105 days at 8-12 km (Table 4). As the altitude increases, the loss 
of 0 3 becomes increasingly controlled by subsidence. Sub- 
sidence from the upper to the lower troposphere over the South 
Atlantic takes place on a timescale of the order of 10 days [Krish- 
namurti et al. , this issue]; we thus expect a characteristic time of 
about 2 weeks for 0 3 to adjust to chemical steady state in the 0 to 
I2-km column. Such a characteristic time is sufficiently long that 
horizontal mass exchange with other regions could represent an 
important term in the budget of tropospheric 0 3 over the South 
Atlantic basin. However, it appears from our close regional bal- 
ance between chemical production and loss that this term does 
not provide either a major net source or a major net sink of 0 3 in 
the South Atlantic basin troposphere, i.e., that the horizontal flux 
divergence is small. A general representation of the 0 3 budget in 
the tropical troposphere consistent with this view is presented in 
section 7. 

Tropospheric 0 3 columns observed over the South Atlantic 
basin during the wet season are about 40% lower than observed in 


TRACE A [Kirchhoff et al. , 1991 ; Olson et al this issue; Thomp- 
son et a! ., this issue]. We examined whether these wet season 
columns could be sustained by in situ photochemical production 
in the same way as in TRACE A. For this purpose we decreased 
the NO and CO levels in the one-dimensional model to represent 
an atmosphere unaffected by biomass burning. Singh et al. 
[1990] reported a median NO concentration of 8 pptv at 0 to 6- 
km altitude over the Amazon basin in April-May (wet season), 
with little vertical gradient. Drummond et al. [1988] reported a 
median NO concentration of 60 pptv at 1 1.5-km altitude between 
Rio de Janeiro and Dakar in June, at the beginning of the dry sea- 
son when biomass burning influence is small and tropospheric 0 3 
columns still assume wet season values [Kirchhoff et al. , 1991], 
We took these observations, which are about 3 times lower than 
in TRACE A, as constraints in the model. We further adopted a 
background CO concentration of 50 ppbv. Hydrocarbon concen- 
trations were not modified from their values in Table la because 
they have little effect on 0 3 production. As shown in Figure 5, 


Table 4. Budget of Tropospheric Ozone Over the South Atlantic Basin 

Altitude, km 


0-4 4-8 8-12 

Number of Model Points 



292 

246 

574 

Production and loss, ppbv d _1 




Pok 

8.0 

3.8 

2.9 

f'Ox 

9.0 

3.6 

0.7 

(P~U Ox 

Column rates, I0 ln cm 2 s -1 

-0.7 

0.3 

2.1 

^Ox 

77 

23 

i 1 


80 

24 

3 

(P-U Ox 

-6 

2 

8 

0 3 lifetime ([0 3 J/L 0x ), days 

6 

25 

105 

0 3 replacement time ([0 3 ]/Po x )» days 

7 

20 

24 


Median 24-hour average rates and characteristic times computed in the model for the southern 
hemisphere TRACE A data. O x is odd oxygen (O x = 0 3 + O + N0 2 + HN0 4 + (2x)N0 3 + 
(3x)N 2 0 5 + peroxyacylnitrates). Column rates refer to the 4-km altitude bands. Because values are 
medians instead of means, (P-L ) ()x and P Gx - L 0x are not identical. 
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0 3 , ppbv 

Figure 5. Steady state concentrations computed in a one- 
dimensional photochemical model constrained with the median 
TRACE A observations for 0 3 precursors (solid line) and with 
NO and CO concentrations reduced to reflect a background atmo- 
sphere free of biomass burning influence (dotted line). The model 
uses a zero-flux boundary condition at 12-km altitude. The T 
symbols are the median 0 3 concentrations in the southern hemi- 
sphere TRACE A data. The "W" symbols are the mean wet sea- 
son (March-May) O^ concentrations observed at Natal Brazil 
[Kirchlioffetal., 1991]. 


the vertical O3 profile computed in the one-dimensional model is 
consistent with the mean wet season ozonesonde observations by 
Kirchhoff et at. [1991] at Natal Brazil. It thus appears that the 
tropospheric 0} column over the South Atlantic basin during the 
wet season can be sustained fully by in situ production, similarly 
to the dry season. Although NO, levels are lower than in the dry 
season, the levels to be sustained are also lower. 

The above analysis implies a strong sensitivity of tropos- 
pheric 0 } over the South Atlantic basin to the supply of NO,. 
We show in Figure 6 the dependence of the net Cb production 
rate on NO, concentrations in the model for points representing 
the median conditions observed in TRACE A at 0 to 4 km, 4 to 8 
km. and 8 to 12-km altitude (Table la). Even at the relatively 
hisih NO, concentrations observed in TRACE A, 0 ? production 
is still NO, limited throughout the 0 to 12 km column. A number 
of authors have noted that NO, levels in the remote troposphere 
are near the turnover point for net O^ production versus net loss 
[Chameides et at., 1987, 1990; Jacob et al , 1992; Liu et at., 
1992; Mauzerall et al. , 1996; Carroll and Thompson , 1995], and 
this is also found in the TRACE A data below 8 km (Figure 6). 
Based on the discussion above, we explain this general result as 
reflecting the tendency of Oj to approach a chemical steady state 
determined by the NO, concentration. 


6. Nitrogen Oxides 

Several authors have examined the primary sources of NO, 
over the South Atlantic basin in TRACE A by correlating the 
observed NO concentrations with chemical variables and with 
back-trajectories [Smyth et at., this issue; Singh et al this issue. 
Talbot et al., this issue; Pickering et al., this issue]. They con- 
clude that biomass burning was a major source of NO, 
throughout the tropospheric column, and that lightning was also a 
major and perhaps dominant source in the upper troposphere. 
Soil emissions would provide an additional source of NO, in the 
continental boundary layer [Jacob and Wofsy , 1988, 1990] but the 
TRACE A data are not well suited for assessing the importance of 
this source. Transport from the stratosphere was a negligible con- 
tributor to the NO, budget even in the upper troposphere [Smyth 
et at., this issue]. 

Chemical cycling between NO, and its oxidation products 
played probably a major role in maintaining high NO, concentra- 
tions over the South Atlantic during TRACE-A. The mean life- 
time of NO, against reactions (8)-(10) computed in the model for 
the TRACE A conditions is only 4 hours at 0-4 km and 2.4 days 
at 8-12 km. Plate 2 shows the distribution of NO, concentrations 
over the South Atlantic basin as a function of longitude and alti- 
tude (NO is measured and NOj is from the photochemical 
model). Concentrations are high in the continental boundary 
layer and the upper troposphere, consistent with primary sources 


Sensitivity of (P-L)03 to NOx, median TRACE-A conditions 


TRACE-A median 




Figure 6. Dependence of the computed 24-hour average net O^ 
production rate on NO, concentrations for the TRACE A condi- 
tions. All variables other than NO, are specified at their median 
TRACE A values given in Table 1 a for 0-4, 4-8, and 8- 1 2 km alti- 
tude. The dashed vertical lines show the median NO, concentra- 
tions in the southern hemisphere TRACE A data for the different 
altitude bands. 
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Plate 2. Mean concentration of NO x in TRACE A at 0°-30°S latitude, as a function of loncitude and altitude. 


from combustion, soils, and lightning; during deep convective 
events over the continents, NO x from combustion and soils is 
injected to the upper troposphere, and additional NO x is produced 
by lightning. However, above 4 km we find no significant 
decrease of NO x from continents to oceans, even though the pri- 
mary sources of NO x are continental (there is little lightning over 
the oceans) and the lifetime of NO x is short. It appears that the 
high concentrations over the oceans must be maintained by chem- 
ical recycling. 

We used our photochemical model to diagnose the degree to 
which NO x levels in TRACE A could be maintained by chemical 
recycling based on current understanding of tropospheric chemis- 
try. For this purpose we computed the chemical production and 
loss rates of NO x at individual points in the merged time series 
matched to the HNO, sampling times, using the model con- 


strained with measurements of NO, HN0 3 , and PAN (details of 
the approach are given in section 2). Figure 7 shows the ratio of 
chemical loss to chemical production of NO x for individual 
points as a function of altitude, and Table 5 summarizes the con- 
tributions of individual reactions to the NO x budget at different 
altitudes. Below 4 km there is considerable scatter for individual 
points, reflecting primarily the relative abundances of NO x and 
PAN; on average, we find a balance between the source of NO x 
from decomposition of PAN and the sinks of NO x from oxidation 
to HNO< and PAN. The data below 4 km are not inconsistent 
with a scenario where PAN produced in biomass burning plumes 
decomposes to sustain NO x levels in the lower troposphere on the 
regional scale, but the paucity of points precludes any reliable 
analysis. 

Above 8 km we find a systematic imbalance in the chemical 
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Figure 7. Ratio of 24-hour average chemical loss to chemical pro- 
duction rates of NO x as a function of altitude for individual meas- 
urement time intervals in the southern hemisphere TRACE A data. 


budget for NO x (Figure 7), with production ot NO, balancing on 
average only 16% of the loss. This deficit is due to insufficient 
recycling of HNO^ to NO x ; recycling of PAN is not an issue 
because rates for NO x -PAN cycling are slow and PAN is near 
chemical steady state (Table 5). The median rate ot conversion 
of NO, to HNO^ at 8-12 km in the model is 57 pptv day' 1 , while 
the median observed HNO^ concentration at that altitude is 59 
pptv. A HNO* lifetime of 1 day would be necessary to balance 
the HNO^ budget, but the median lifetime of HN0 3 against pho- 


tolysis and reaction with OH is 12 days. It thus appears that a 
major mechanism for high-altitude conversion of HNO, to NO, 
may be missing from the model. The problem would remain 
even if N 2 0^ conversion to HNO^ in aerosols did not proceed, as 
this reaction accounts, on average, for only 30% of the conver- 
sion of NO, toHNO., (Table 5). 

The above analysis suggests that HNO-, in the upper tropo- 
sphere must be consumed on a timescale of 1 day to eventually 
regenerate NO x . Chatfield (1994) proposed that rapid reaction of 
HNO} with CH 2 0 in acid aerosols might take place and yield 
NO, and formic acid as immediate products. We find, however, 
no correlation between HNO^ and NO concentrations in the 
TRACE A data above 8 km (Figure 8), which argues against 
direct conversion of HNO^ to NO, on a short time scale, han et 
al. [1994] proposed that the reaction of HNO-, with CH : 0 in acid 
aerosols may not produce NO x but hydroxymethylnitrate 
H,C(0H)0N0 2 ; this species could then photolyze to yield NO, 
on a timescale of one week, based on laboratory data for other 
organic nitrates [Barnes et al ., 1993]. The Fan et al. [1994] 
mechanism would be more consistent with the lack ot correlation 
between NO and HNO^ in the TRACE A data. 

7. Schematic for Ozone in the Tropical Tropo- 
sphere 

Our analysis of the TRACE A data suggests a simple 
schematic for the regional budget of tropospheric O^ over the 
South Atlantic basin. Continental NO x emissions from South 
America and southern Africa are pumped to the middle and upper 
troposphere by convection, leading to ubiquitous net production 
of o, in the upper troposphere. Eventually, the air subsides over 
the Atlantic Ocean and net loss ot 0-* takes place. Low-level 
adveclion to the continents closes the 0-* cycle. We propose that 
this schematic can be extended to other tropical regions and pro- 
vides a general explanation for the origin of O^ in the tropical tro- 
posphere, as illustrated in Figure 9. In this figure, the large-scale 
circulation of the tropical atmosphere is represented by the 


Table 5. Chemical Cycling of NO x over the South Atlantic Basin 




Altitude Band, km 



0-4 

4-8 

Number of Model Points 

8-12 


15 

37 

149 

NO x loss, pptv d _l 
NOi + OH 

230 

50 

41 

N->Os + aerosol 

141 

37 

16 

CH,CO, + NO, 

132 

19 

7 

RO-i + NO, a 

74 

4 

1 

67 

Total Loss 

787 

112 

NO x production, pptv d " 1 
HNO} photolysis 

14 

3 

i 

HNO 3 + OH 

24 

3 

- 

PAN thermolysis 

728 

6 

<1 

PAN photolysis 

6 

8 

7 

Total Production 

900 

3 1 

1 1 

NO x lifetime, days 

0.18 

0.66 

2.4 


Median 24-hour average rates and lifetimes computed in the model for the southern hemisphere 
TRACE A data. Because the values shown are medians, the sum of the individual rates does not add 
up to the total production or loss of NO x , and the median lifetime does not match exactly the 
corresponding concentration and loss rate; this effect is particularly noticeable at 0-4 km where the 
scatter of rates for individual points is large and the number of points is small. The rates ot produc- 
tion and loss of PAN have been corrected for internal cycling within the (CH^CO.i + PAN) family as 

described in Appendix B. . 

a Formation of all organic nitrates other than PAN. These organic nitrates are assumed men m the 

model. 
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Figure 8. Relationship of HN0 3 and NO concentrations in the 
TRACE A daytime data (0800-1600 local time) at 8 to 12-km alti- 
tude. 


Our representation of the budget of 0 3 in the tropical tropo- 
sphere goes back to the earliest papers arguing that the tropos- 
pheric 0 3 budget is dominated by in situ production and that tran- 
sport from the stratosphere is small in comparison [Chameides 
and Walker , 1973; Crutzen , 1973; Fishman et al, 1979; Liu etal 
19801. Much uncertainty was attached to these early studies 
because of the lack of reliable NO data. Liu et al. [1980] first 
proposed the idea of a photochemical steady state of 0 3 in the 
tropospheric column with net production in the upper troposphere 
compensating for net loss in the lower troposphere. They argued 
that the NO* in the upper troposphere driving the in situ produc- 
tion would originate from the stratosphere, but the TRACE A 
data as well as previous work [Kasibhatla et al ., 1991; Penner et 
al., 1991; Murphy et al. , 1993] show that stratospheric influx is 
negligible compared to combustion and lightning as a source of 
NO x in the tropical upper troposphere. 

8. Conclusions 


Walker cells [after Newell , 1979], Transfer of air from the lower 
to the upper troposphere takes place mainly by deep wet convec- 
tion over South America, Africa, and Oceania (rising branches of 
the Walker cells). These convective motions inject to the upper 
troposphere NO x emitted over the continents by combustion and 
soils, and additional NO* produced by lightning in the convective 
column. The injection of NO x then drives 0 3 production in the 
upper troposphere. It must be pointed out that although produc- 
tion of 0 3 in the upper troposphere is necessary to explain the 
observed increase of 0 3 concentrations with altitude and to bal- 
ance the net 0 3 loss in the lower troposphere, most of the 0 3 pro- 
duction in the tropospheric column actually takes place in the 
lower and middle troposphere (Figure 4). 

Some support for the schematic in Figure 9 is offered by data 
from the equatorial Pacific. Piotrowicz et al. [1991] found that 
the distribution of 0 3 concentrations in surface air over the equa- 
torial Pacific is highly sensitive to seasonal and interannual move- 
ments of the Pacific Walker cell. Fishman et al. [1990] reported 
a west-to-east gradient of increasing tropospheric 0 3 columns 
across the Pacific which may be explained by net production 
in the upper troposphere westerlies and net loss in the lower tro- 
posphere easterlies. It is well established that the lower tropo- 
sphere of the equatorial Pacific is a net sink for 0 3 [Liu et al., 
1983; Johnson et al., 1990]. Limited aircraft data for the western 
tropica] Pacific indicate net 0 3 loss taking place below 6-km alti- 
tude and net production taking place above, with production 
balancing 60% of 0 3 loss in the 0 to 12-km column [Davis et al., 
1996]. 


Photochemical modeling of TRACE A aircraft observations 
over the South Atlantic basin up to 12-km altitude indicates a 
close balance between chemical production and loss of 0 } in the 
0 to 12-km column, with net production in the upper troposphere 
balancing weak net loss in the lower troposphere. Influx from the 
stratosphere is unimportant as a source of tropospheric 0 3 in the 
region. We conclude that the tropospheric 0 3 column can be 
maintained solely by in situ photochemistry. The photochemical 
production of 0 3 is limited by the supply of NO x , even in the 
upper troposphere. Past observations of the remote troposphere 
indicate that NO x concentrations are generally near the turnover 
point for net 0-, production versus net 0 3 loss; we interpret this 
observation as reflecting the tendency of 0 3 concentrations to 
adjust to a chemical steady state defined by the available NO x . 

Biomass burning and lightning were the two primary sources 
of NO x over the South Atlantic basin in TRACE A. Considering 
the short lifetime of NO x against oxidation (only 2-3 days at 8-12 
km altitude), it appears that the high NO* levels observed 
throughout the South Atlantic basin were maintained by chemical 
recycling of NO* from its oxidation products. The low 
HN0 3 /NO concentration ratios observed above 8 km imply rapid 
conversion of HN0 3 to NO x by a mechanism not considered in 
current models. 

We presented a simple schematic for the origin of 0 3 in the 
South Atlantic basin troposphere linking the large-scale Walker 
circulation to in situ photochemistry. Convection over South 
America and southern Africa injects NO x from combustion, soils, 
and lightning to the middle and upper troposphere, leading to net 



Figure 9. Schematic depiction of the O, budget in the tropical troposphere. Net O, production in the upper tropo- 
sphere balances net O, loss in the lower troposphere. Mass exchange between the lower and the upper troposphere is 
represented by the Walker circulation [after Newell , 1979]. Deep convective motions over South America, southern 
Africa, and Oceania (rising branches of the Walker circulation) inject NO, from combustion, soils, and lightning to the 
upper troposphere, driving O, production. Eventually, the air subsides over the oceans and net O, loss takes place in 
the lower troposphere due to low NO x concentrations and high humidities, closing the O, cycle. 
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O; production; eventually, the air subsides over the Atlantic 
Ocean and net loss of O, takes place in the lower troposphere due 
to reduced NO x levels and high humidities. Low-level advcction 
of marine air to the continents closes the 0 3 cycle. We proposed 
that the same schematic applies to other tropical regions and 
offers a general representation of the budget in the tropical 
troposphere. This mechanism would imply a great sensitivity ot 
the global oxidizing power of the atmosphere to NO x emissions 
from tropical continents. 

Peroxide concentrations measured aboard the aircratt in 
TRACE A offered a sensitive test of the CL production and loss 
rates computed in our photochemical model. We found that the 
model captures 78% and 74% of the observed variances of H 2 0 2 
and CPLOOH concentrations, respectively. Below 8-km altitude 
there is no mean bias for H 2 0 2 and a 30% negative bias for 
CFLOOH. Above 8-km the model underestimates H 2 0 2 by a 
factor of 1 .8 and CH,OOH by a factor of 3.6. Even allowing for 
these discrepancies, the agreement between modeled and 
observed peroxide concentrations is sufficiently good to lend 
confidence in our conclusions regarding the CL budget in the tro- 
pospheric column. 

Formaldehyde concentrations measured aboard the aircraft 
provided another test of photochemistry in the model. The model 
captures 62% of the variance of measured CH 2 0 concentrations 
but tends to overestimate concentrations in the lower troposphere 
(0-4 km) and underestimate concentrations in the middle tropo- 
sphere (4-8 km). The largest discrepancy between model and 
observations is in the marine boundary layer, where the model is 
too high by a factor of 3. A review of previous data in the tropi- 
cal marine boundary layer indicates considerable variability in 
measured CH.O concentrations that cannot be explained with 
current models. The issue requires further investigation. 

Appendix A: Photochemical Model 

Our photochemical mechanism includes inorganic and 
methane chemistry from DeMore et al. [1994], and non-methane 
hydrocarbon chemistry from Atkinson et al. [1993] and Atkinson 
[1994]. The rate constant and branching ratio for the 

CH^OOH+OH reaction are from Vcighjiani and Ravishankara 
[1989]. The yield of O('D) from photolysis of 0 4 is from 
Michelsen et al. [1994], Recent data arc used for the photolysis 
of HNO, [Burkholder et al., 1993], PAN [Talukdar et al., 1995]. 
and ketones [Raher and Moortgaat , 1995]. 

Hydrolysis of N 2 0 5 to HNCL in aerosols is included with a 
rate constant determined by the collision frequency ot N 2 0 5 
molecules on aerosol particles [Fuchs and Sutugin , 1971] and a 
reaction probability of 0.1 based on data for aqueous surfaces 
| DeMore et a!., 1994]. The collision frequency is computed 
using the aerosol size distribution (range 0.12 to 2 pm radius) 
measured aboard the aircraft with 10-s resolution. The reaction 
probability of N 2 0 5 on the aerosol surface would be low if the 
aerosol were solid [Mozurkewich and Calvert , 1988). At the tem- 
peratures and humidities found above 8-km altitude in TRACE A. 
the stable phase of the binary H 2 S0 4 -H 2 0 system is usually solid 
H 2 S0 4 *4H 2 0 (sulfuric acid tetrahydrate, or SAT). The reaction 
probability of N 2 CL on SAT is only 0.001-0.01 [Hctnson and 
Ravishankara. 1993]. However, the aerosol could be in a meta- 
stable liquid state. Considering that our model uses NO as an 
input variable, the inclusion ot N 2 0^ hydrolysis has no significant 
effect on the model photochemistry, but it provides a significant 
term in the computed mass balance for NO x . On the one hand, 
the rale of hydrolysis could be underestimated it there is 
significant aerosol surface area outside of the 0.12 to 2-pm parti- 
cle size range; on the other hand, it could be overestimated if the 
aerosol includes a significant population of dry particles. 

Photolysis frequencies are computed at each point using the 
clear-sky radiation model of Logan et al. [ 1 98 1 ] with a surface 


albedo of 0.1. They are then scaled for cloud and albedo effects 
using a factor CF = 7 N()2 . F pp k . y /,/ N o 2 . t i ca . , where J^n 2X \ Cdr is the 
clear-sky NO : photolysis frequency computed with the radiation 
model at the local time of observation, and ./ N( , 2 . h pp i^ is the value 
derived from the Eppley UV radiometer measurements using the 
parameterization of Madronich [1987], as described by Davis et 
al. | 1993]. This parameterization is applicable only to solar zen- 
uh angles less than 70° and we therefore restrict our calculations 
to this range (97% of observations). Mean CF values in TRACE 
A were 0.99 at 0-2 km and 1.35 at 8-12 km; values of CF higher 
than 2 (maximum 2.51) were found when the aircratt fiew 
through clouds or above stratus decks under high-Sun conditions. 
There was no significant dependence of CF on Sun angle. To 
check for consistency between the va ^ ues computed from 
our clear-sky radiation model and those inferred from the Eppley 
data, we searched for cloud-free conditions in the aerosol lidar 
data and in videos taken from forward and nadir cameras aboard 
the aircraft. The only unambiguous clear-sky period was on flight 
12 over Namibia (1200-1240 GMT), with the aircraft cruising at 
10-km altitude. The mean CF for this period was 1.06. 

Model results for CH^OOH and CH 2 0 are sensitive to the 
oxidation mechanism assumed for acetic acid. Acetic acid con- 
centrations measured in TRACE A averaged 2-3 ppbv throughout 
the 0 to 12-km column (Table la). No laboratory data are avail- 
able to our know ledge for the products of acetic acid oxidation by 
OH. The mechanism proposed by Atkinson [1994], 

O, 

CHOCOOH + OH — >CH,0 2 +C0 2 + H : 0 (Al) 


would result in acetic acid providing a major source of CH^OOH 
and CH 2 0 in the upper troposphere; at 230 K, 3 ppbv of acetic 
acid would produce as much CH 2 0 as 1700 ppbv of CH 4 . How- 
ever, the TRACE A data show no correlation between measured 
acetic acid and CH.O in the upper troposphere, even when meas- 
ured acetic acid is as high as 8 ppbv (Figure A I ). We conclude 
that either there is an error in the measurements, or the oxidation 
of acetic acid does not produce CH 2 0. To side step this problem, 
we assume for the purposes of the model that the oxidation of 
acetic acid follows the scheme of Madronich and Calvert [ 1989], 
which bypasses CH^OOH and CH 2 0 formation: 

CHOCOOH + OH — > CH 2 (0 2 )C00H + H 2 0 (A2) 

CH 2 (0 2 )COOH + NO — > CH 2 (0)COOH + NO : (A3) 

CH 2 (0)C00H + 0 2 CHOCOOH + H0 2 (A4) 

CHOCOOH + OH — > + CO + C0 2 + H.O + HO : ( A5) 


0.6 
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, n ’ ' • 


, iy * 

0.0 



0 2 4 6 8 

ACETIC ACIO, ppbv 

Figure Al. Relationship of CH 2 0 and acetic acid concentrations 
in the TRACE A data at 8 to 12-km altitude. 
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Figure A2. Relationship of acetic acid and CO concentrations in 
the southern hemisphere TRACE A data (squares, crosses) and in 
the DECAFE aircraft data over the Congo rainforest [Helas et aL , 
1992]. The solid and dashed lines are linear best fits to the 
TRACE A and DECAFE data, respectively. The DECAFE data 
set consists of only 10 points. 


(13) and (14) to maintain the CH 3 C0 3 steady state. Use of 
uncorrected rates R m , R u , and R {4 to construct a NO x budget can 
lead to the illusion of a ciose balance between production and loss 
of NO x and to an underestimate of the chemical lifetime of NO x . 

We account for this effect in our NO x budget analyses by 
using corrected rates R'\ 


tf'n, = -Kin 
P'\y - R\j 

*'l4=*14 


P PANx + + R \ 4 

^PANx 

^PANx + ^13 + R\4 

^PANx 

P PANx + P\y + ^ !4 


(Bl) 

(B2) 

(B3) 


where an ^ ^panx are the chemical production and loss rates 
of the PAN X family. In the TRACE A data set, the correction 
factor PpanxAPpanx + P\y + P \4 ) ranges from 0.03 to 0.72 at 0-4 
km and from 0.80 to 0.99 at 8-12 km. The correction factor 
^PANx^PpANv + P\y + R ranges from 0.23 to 0.71 at 0-4 km and 
from 0.76 to 0.97 at 8-12 km. The largest corrections are in the 
lower troposphere, where thermal decomposition of PAN often 
dominates the CH^COj source. 


When NO concentrations are low, formation of the organic 
peroxide CH 2 (OOH)COOH provides a temporary reservoir but 
does not change the ultimate fate of the carbon. With this 
mechanism, acetic acid has little effect on the model photochem- 
istry even at the high concentrations measured in TRACE A. 

The origin of the high acetic acid concentrations observed in 
TRACE A remains a puzzle. There are to our knowledge no pre- 
vious observations in the upper troposphere that can be used for 
comparison. Figure A2 shows the relationship between acetic 
acid and CO in TRACE A, and compares to the relationship 
found in 10 dry season measurements taken at 0-4 km over the 
Congo rainforest during the Dynamique et Chimie en Foret Equa- 
torial (DECAFE) expedition [Helas et aL, 1992]. Both the 
TRACE A and the DECAFE data show correlations between 
acetic acid and CO at 0-4 km, with similar slopes, indicating 
acetic acid emission from biomass burning; however, the TRACE 
A data include, also, a major component not correlated with CO, 
resulting in acetic acid concentrations about 3 ppbv higher than in 
DECAFE for the same amount of CO. The acetic acid concentra- 
tions measured in the upper troposphere in TRACE A were not 
correlated with CO or any of the other variables measured aboard 
the aircraft except formic acid. They do, however, show a pattern 
of higher values over Brazil and Africa than over the South 
Atlantic. 

Appendix B: Treatment of NO x -PAN Inter- 

conversion in NO x Budget Calculations 

When assessing the importance of NO x -PAN interconversion 
in a regional NO x budget calculation constrained with observa- 
tions for NO x (or NO) and PAN, and with CH^CO-* taken to be at 
steady state, one must correct the rates of reactions (10), (13), and 
(14) to account for cycling within the chemical family PAN X = 
CHjCOj + PAN. To illustrate the need for this correction, con- 
sider a situation where reactions (10), (13), and (14) dominate the 
production and loss of CH^CO*. Under this situation, given any 
specified input pair of NO x and PAN concentrations, the steady 
state CHjCO* concentration computed in the model yields an 
apparent steady state between NO x and PAN; this steady state 
plays, however, no role in regulating NO x since loss of NO x and 
CH3CO3 by reaction (10) forces production of NO x by reactions 
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